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We have performed self-consistent Skyrme Hartree-Fock-Bogolyubov calcula- 
tions for nuclei close to ^^^No. Self-consistent deformations, including /32,4,6,8 
as functions of the rotational frequency, were determined for even-even nuclei 
246,248,250pjjj_ 252,254^^^^ 256j^f rpj^^ quasiparticle spectra for N = 151 

isotones and Z = 99 isotopes were calculated and compared with experimental 
data and the results of Woods-Saxon calculations. We found that our calcu- 
lations give high-order deformations similar to those obtained for the Woods- 
Saxon potential, and that the experimental quasiparticle energies are reason- 
ably well reproduced. 
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1. Introduction 

Nuclei in the nobehum region represent the heaviest systems for which de- 
tailed spectroscopic information on the structure is available.^ Hence, they 
provide an excellent testing ground for various theoretical models that aim 
at reliable predictions of properties of superheavy elements. Although signif- 
icant progress has been made both experimentally and theoretically, there 
are still challenges for self-consistent mean-field calculations (see Ref. 2 for 
a recent review). One of the challenges is the correct description of de- 
formed shell gaps. Macroscopic-microscopic calculations predict shell gaps 
at = 152 and Z = 100, whereas typical self-consistent mean-field cal- 
culations give shell gaps at = 150, and Z = 98 and Z = 104. The ex- 
perimental findings are consistent with the macroscopic-microscopic shell 
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gaps. 

Current advanced readjustements"^ of coupling constants of the Skyrme 
functional do not cure this deficiency, although they point to a large under- 
determination of the spin-orbit properties. To study detailed spectroscopic 
properties of heavy and superheavy nuclei, one thus may attempt a fine 
tuning of the spin-orbit coupling constants, so as to bring the deformed 
shell gaps in closer agreement with data. For this purpose, we recently 
readjusted^ the spin-orbit coupling constants, C^'^ and CY'\ and pairing 
strengths, Vq and Vq , of the Skyrme UNEDFl functional^ to match the 
experimental excitation energies of the ll/2~ state in ^^^Cf and 7/2+ state 
in ^^^Bk. For the time-odd coupling constants of the functional, we adopted 
the prescription based on the Landau parameters, as defined in Rcf. 5. The 
remaining parameters of UNEDFl were kept unchanged. We call this new 
parameter set UNEDFl Details of the adjustment procedures will be 
presented in a forth-coming publication.* Here we only quote the obtained 
values, which read, 

{C^'' , C7'0 = (-88.05040, 8.45838) MeVtm^, (1) 

{K^ V) = (-205.05, -252.48) MeV. (2) 

As a result, in Ref. 4 we found improvements of shell gaps, quasipar- 
ticle energies, and moments of inertia (Mol) in nuclei around ^^*No. Also 
a good agreement with experiment has been obtained for the quasiparti- 
cle energies in ^"^^Bk and ^^^Cf. The single-particle spectra more closely 
resemble those obtained with the WS potential, and shell gaps are opened 
at Z = 100 and N = 152 in ^^*No. In the present work, we extended the 
calculations to other miclei in the nobeliiim region. We investigated the 
self-consistent deformations as a function of frequency and we performed 
blocked calculations for odd-j4 N = 151 isotones and Es {Z = 99) isotopes. 

2. The model 

In this work, all the calculations were performed by using the symmetry- 
unrestricted solver HFODD^ (v2.52j). For rotational calculations, the stan- 
dard cranking term —coJy was added to the Hamiltonian. The effect of an 
odd nucleon was treated by employing the blocking approximation, that is, 
by using the density matrix and pairing tensor in the form,^ 

P'fnn = {V*V^)mn + Ur^^U^^ - ^l^K^, (3) 
^mn ~ (^ ^ )mn + UmiJ,Vnn ~ Kn/i^"/^- i^) 

We refer the reader to Ref. 8 for details. 
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In the cranking-plus-blocking calculations, the timc-rcversal symmetry 
is broken, so the time-odd terms of the Skyrme density functional are 
present. The importance of time-odd terms in describing collective rota- 
tional bands has been emphasized in Ref. 9. The advantage of using a 
symmetry-unrestricted solver is that the effects of non-zero time-odd po- 
tentials can be taken into account in a strict manner. 

We used 680 deformed harmonic-oscillator basis states, with the 
basis-deformation parameters of huix = hojy = 8.5049 McV and Hojz = 
6.4823 MeV. The cutoff energy in the quasiparticle spectrum was Ec = 
60MeV. In the pairing channel, to describe the approximate particle- 
number projection, we adopted the Lipkin-Nogami formalism. 

3. Results 

In the nobelium region, high-order deformations of the WS potential were 
found to play important roles in explaining many high-spin phenomena. 
For example, the decrease of i?a;(2+) in ^^^No^° and the delayed upbending 
of Mol in ^^''No as compared to ^^^No^^ were attributed to the effect of /Je 
deformation on the relevant single-particle levels. 

We performed cranking calculations for the 6 nuclei, in which the 
ground-state rotational bands arc observed, namely, for 246,248,250pj-|^^ 
252,254]vjq^ and ^^^Rf. We determined the Mols, proton and neutron pair- 
ing gaps, and self-consistent deformations as functions of the rotational 
frequency.^ The experimental kinematic and dynamic Mols were well re- 
produced by our calculations. Deformations obtained within our fully self- 
consistent calculations provide a reasonable description of polarization ef- 
fects exerted by the occupations of individual deformed orbitals. They were 
determined as the Bohr deformation parameters Px of sharp-edge mass dis- 
tributions that have multipole moments Qx, for A = 2-8, identical to those 
calculated microscopically, see Ref. 12. All nuclei considered here have axial 
ground-state shapes. 

Figure 1 shows deformations /32.4,6,8 as a function of rotational frequency 
w. At w = 0, deformations (32 are systematically higher than those obtained 
by macroscopic-microscopic WS calculations, and, in general, lower 
than those obtained within other rcilativistic and non-relativistic mean-field 
calculations.^'''^'' It is interesting to note that the high-order deformations 
presented in Figure 1 are very close to those of Refs. 10 and 11. With in- 
creasing rotational frequency, the deformations stay almost constant. Only 
at uj ^ 0.3 MeV, deformations 02 begin to decrease. This behavior is con- 
sistent with that of Ref. 11. 
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Fig. 1. Deformations /92,4,6,8 a-s functions of the rotational frequency, calculated self- 
consistently for six nuclei,' 246,248, 250pjn, 252,254^^^ 256j^f 



Figure 2 (left) shows the calculated quasiparticle spectra for N = 151 
isotones from Z = 94 to 104. For these nuclei, experimental data are avail- 
able in this mass region. The observed ground states are dominated by 
the 9/2[734] orbital, with the 5/2[622] state located at ~200keV. In ^^Spu 
and ^'*'^Cm, our calculations predict ground states of 9/2[734]. For heavier 
isotones, the 7/2 [624] state becomes lower than the 9/2 [734] state, but only 
by several tens of keV. These two levels are right below the N = 152 shell 
gap. They are very close to each other, and change order with increasing 
proton number. This indicates that (i) there is a significant TV-dependence 
of the single-particle and quasiparticle energies and (ii) adjusting the spin- 
orbit force to obtain better agreement for the ll/2~ state in ^^'^Cf results 
in the ji5/2 and gg/2 spherical shells being too close together. The present 
calculations predict a 5/2 [622] state which is too high and underestimate 
the energies of the 7/2[624] and l/2[620] states. 

Comparing our results with those of Ref. 16, we note the presence of 
low-lying 7/2 [624] states there too. We also note an overestimation of their 
1/2 [620] energies. The excitation energy of the 1/2 [620] orbital is related 
to the shell gap at = 152. Its isospin dependence is very interesting - we 
note that our results reproduce the observed increase of excitation energy 
of 1/2 [620] with increasing proton number. 

Experimental information on quasi-proton spectra are relatively scarce.^ 
Close to Z = 100, systematic spectra are only observed for Es {Z = 99) 
isotopes. Figure 2 (right) displays the calculated quasiparticle spectra for 
Z = 99 isotopes from N = 144 to 154. One can see that for the ground states 
the calculations predict two very close states, 3/2[521] and 7/2[633], which is 
consistent with experimental observations.^^ Compared to the calculations 
presented in Ref. 17, our results show a much lower 3/2[521] orbital. This 
is because our calculations predict very close 7/2[633] and 3/2[521] levels 
below the Z — 100 shell opening, whereas the WS calculations give them 
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Fig. 2. Quasiparticle energies calculated for the UNEDFl|^*^ parameter set,* and for 
the N = 151 isotones (left) and Es [Z = 99) isotopes (right). 

~200keV apart. The 7/2 [514] state is also systematically observed in the 
Es isotopes. However, our calculation fails to reproduce the increase of 
excitation energy of this orbital with increasing neutron number. We also 
note that in our calculations, the 5/2 [523] states exhibit an increase similar 
to that observed in experiment for the state attributed to the 7/2 [514] 
orbital. Further experimental and theoretical work is needed for these odd- 
Z nuclei. 

4. Conclusion 

We performed self-consistent Skyrme Hartree-Fock-Bogolyubov calcula- 
tions for nuclei in nobelium region. Deformations /32,4,6,8 were studied for 
even-even nuclei 246,248,250pjjj^ 252,254]-^^^ 256j^f_ r^j^g quasiparticle spec- 
tra for N — 151 isotones and Z = 99 Es isotopes were systematically 
calculated by using the blocking approximation within the Hartree-Fock- 
Bogolyubov method. Comparisons were made with the available experi- 
mental data and WS calculations. It was found that our results give similar 
deformations to those obtained within the WS calculations. The observed 
quasi-particle energies are reasonably well reproduced, especially for proton 
quasiparticle states. 
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